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Particle production and event topology are very strongly correlated in high-energy hadronic and 
nuclear collisions. Event topology is decided by the underlying particle production dynamics and 
medium effects. Transverse spherocity is an event shape observable, which has been used in pp and 
heavy-ion collisions to separate the events based on their geometrical shapes. It has the unique 
capability to distinguish between jetty and isotropic events. In this work, we have implemented 
transverse spherocity in Pb-Pb collisions at ./snn = 5.02 TeV using A Multi-Phase Transport Model 
(AMPT). While awaiting for experimental explorations, we perform a feasibility study of transverse 
spherocity dependence of some of the global observables in heavy-ion collisions at the Large Hadron 
Collider energies. These global observables include the Bjorken energy density (¢)), squared speed of 
sound (c2) in the medium and the kinetic freeze-out properties for different collision centralities. The 
present study reveals about the usefulness of event topology dependent measurements in heavy-ion 
collisions. 


Heavy-ion collisions at the ultra-relativistic energies aim to produce a deconfined state of quarks and gluons, the 
primordial matter believed to have formed at the infancy of the Universe. The matter created in such collisions at 
the Relativistic Heavy-Ion Collider (RHIC) at the Brookhaven National Laboratory, USA and at the Large Hadron 
Collider (LHC) at European Center for Nuclear Research (CERN), Switzerland gives an opportunity to study its 
properties at the extreme conditions of temperature and energy densities. Global properties of the created matter 
such as total charged particle multiplicity in the final state, initial energy density, and temperature of the system 
play a pivotal role to understand the form of the created matter, while addressing many fundamental questions 
in basic science. The expansion of the created fireball because of huge concentration of initial energy density 
and high temperature is probed through the equation of state and hence the speed of sound in the medium. The 
final state particle abundances are governed by the initial state of the matter—if the collision creates a partonic 
deconfined colored phase of matter or a confined, color neutral hadronic state. The event topology is governed 
by the underlying particle production mechanism. For example, a back-to-back momentum conserving shower 
of particles, called jets is an event topology, whose underlying mechanism is governed by hard perturbative 
Quantum Chromodynamics (pQCD) high transverse momentum processes, whereas an isotropic geometry in 
the event topology is mostly rich in soft QCD non-perturbative interactions. 

In elementary and hadronic collisions at GeV and TeV energies, although event topology dependent stud- 
ies have got some level of importance, because of a dense medium formation in heavy-ion collisions, these 
techniques like sphericity, transverse spherocity, Rr etc. are not applied to heavy-ion collisions. Event topology 
dependent characterisation of the systems produced in heavy-ion collisions through the global properties are 
not a well explored area. This is a first attempt to explore the sensitiveness of the global observables in heavy-ion 
collisions, to the event topology and hence the underlying particle production dynamics. In our recent study', we 
found that the anisotropic flow strongly depends on transverse spherocity in heavy-ion collision systems. Thus, 
it would be interesting to see how the global properties, often studied in heavy-ion collisions, vary as a function 
of transverse spherocity. It is worth to note that global properties like Bjorken energy density and speed of sound 
give insights to the initial state of the produced system while the kinetic freezeout parameters provide insights to 
the evolution of produced particles in the medium. The use of transverse spherocity along with collision central- 
ity also provides an opportunity to study such observables in a multi-differential way. In the present work, we 
use transverse spherocity as the event shape observable and study the global properties such as Bjorken energy 
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density, speed of sound and kinetic freeze-out parameters for different centrality classes in Pb-Pb collisions at 
./SNN = 5.02 TeV using A Multi-Phase Transport Model (AMPT). 

‘The paper is organised as follows. We begin with a brief introduction about the usefulness of event topology 
studies in heavy-ion collisions. In “Event generation and analysis methodology’, the event generation methodol- 
ogy in AMPT and the definition of transverse spherocity is given. We report and discuss the results in “Results 
and discussion”. The results are summarised in “Summary and conclusion”. 


Event generation and analysis methodology 
In this section, we begin with a brief introduction on AMPT model. Then, we proceed to define the transverse 
spherocity as an event shape analysis tool. 


A Multi-phase Transport (AMPT) Model. A Multi-Phase Transport Model contains four components 
namely*6°"°, 


e Initialisation of collisions using HIJING model: the cross-section of the produced mini-jets in pp collisions 
is calculated and then converted to heavy-ion collisions via inbuilt Glauber model 

e Parton transport after initialisation: transportation of produced partons is performed via Zhang’s parton 
cascade model 

e Hadronisation mechanism: in string melting version, the transported partons are hadronised using spatial 
coalescence mechanism; in the default AMPT version, fragmentation mechanism using Lund fragmentation 
parameters are used for hadronising the transported partons 

e Hadron transport: the hadrons undergo evolution in relativistic transport mechanism via meson-baryon, 
meson-meson and baryon-baryon interactions 


As, the particle flow and spectra at the mid- pr regions are well explained by quark coalescence mechanism for 
hadronisation®””°, we have used string melting mode for all of our calculations. We have used the AMPT ver- 
sion 2.26t7 (released: 28/10/2016) in our current work. The AMPT settings in the current work, are the same as 
reported in Ref.'""!. For the input of impact parameter values for different centralities in Pb-Pb collisions, we have 
used Ref.!*. One should note here that, high centrality collisions correspond to low impact parameter values and 
higher final state charged-particle multiplicity ((dNc,/dn)). Although the concept of centrality is widely used in 
heavy-ion collisions, in view of a final state multiplicity scaling across collisions species, that is observed at the 
LHC energies, we may use centrality and (dN.4,/dn) variably in this work. 


Transverse spherocity. Transverse spherocity is defined for a unit vector f(ny,0) that minimizes the 


ratio!>"; 

n? (Si |pr, x fl \* 

So = —( =}. (1) 
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fis an arbitrary unit vector in the transverse plane. To find such a unit vector, one has to perform iteration 
through all possible values of # with azimuthal angle 0 to 27 and select the # such that the term inside the bracket 
given in Eq. (1) becomes minimum for a given event. Here, the index i runs over all the final state hadrons in an 
event. By construction, transverse spherocity is infrared and collinear safe!’ and the extreme limits are related 
to specific configurations of events in transverse plane. Transverse spherocity becoming 0 are the events with 
pencil-like (back-to-back) structure and called as jetty events while 1 would mean the events are isotropic. A 
schematic picture illustrating the event topology is shown in Fig. 1. For the sake of simplicity, here onwards, the 
transverse spherocity is referred as spherocity. The spherocity distributions are selected in the pseudorapidity 
range of|n| < 0.8 with a minimum constraint of 5 charged particles with py > 0.15 GeV/c to recreate the similar 
conditions as in ALICE experiment at the LHC. The jetty events are those events having spherocity values in the 
lowest 20 percent and the isotropic events are those occupying the highest 20 percent in the spherocity distribu- 
tion of all the events. The spherocity cuts for each centrality are mentioned in Table 1. 

Top panel of Fig. 2 shows the spherocity distributions for different centrality classes in Pb-Pb collisions, 
/SNN = 5.02 TeV at mid-rapidity (/n| < 0.8) using AMPT model and similarly, the bottom panel is the spheroc- 
ity distribution for different charged particle multiplicity classes in pp collisions, ,/s = 13 TeV at mid-rapidity 
(In| < 0.8) using PYTHIA8"*. The details of event generation methodology using PYTHIA8 can be found in 
Ref.'”. For pp collisions, the charged particle multiplicities are chosen in the acceptance of V0 detector in ALICE 
at the LHC with pseudorapidity coverage of VOA (2.8 < n < 5.1) and VOC (—3.7 < n < —1.7). One should 
note here that the particle production mechanisms in AMPT and PYTHIA8 models are completely different 
but we have chosen the tunes of the models where the models describe many of the experimental observables. 
At a first glance, the spherocity distributions in Pb-Pb collisions look shifted more towards the isotropic limit 
when compared to the pp collisions, where the distributions are shifted towards the jetty limit. This behavior is 
understood based on the fact that the system size in Pb-Pb collisions are significantly higher when compared 
with pp collisions and the medium effect in terms of rescattering helps taking the system towards isotropisa- 
tion. In comparison to pp collisions, where we observe a distribution of spherocity or in other words, there is 
an equal production probability of both jetty and isotropic events, medium effects in heavy-ion collisions in 
principle destroy the jettiness event topology, which is seen from Fig. 2. A skewed distribution of spherocity 
towards isotropic limit is an indication of the formation of a QCD medium in heavy-ion collisions. When stud- 
ied as a function of centrality (multiplicity) classes for Pb-Pb (pp) collisions, it is observed that the spherocity 
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Isotropic (Sy > 1) 


Figure 1. Schematic picture showing jetty and isotropic events in the transverse plane, assuming the z-axis is 
the beam axis or the longitudinal axis. 





























Centrality (%) | Low-So | High-So 
0-10 0-0.880 | 0.953-1 
10-20 0-0.813 | 0.914-1 
20-30 0-0.760 | 0.882-1 
30-40 0-0.735 | 0.869-1 
40-50 0-0.716 | 0.865-1 
50-60 0-0.710 | 0.870-1 
60-70 0-0.707 | 0.873-1 
70-100 0-0.535 | 0.822-1 

















Table 1. Low 20 % and high 20% cuts on spherocity distribution in Pb-Pb collisions at ./syn = 5.02 TeV for 
different centrality classes. 


distributions are shifted towards the isotropic limit for central (high-multiplicity) collisions compared to periph- 
eral (low-multiplicity) collisions. 

Now, we proceed to discuss the global properties and their dependence on spherocity classes in the next sec- 
tion. For the sake of simplicity, here onwards we refer 7+ + 2~,K*+K , and p+pas pions, kaons, and protons, 
respectively. 


Results and discussion 

Bjorken energy density (€gj). In heavy-ion collisions, the transverse energy (Er) is one of the significant 
global observables that is used to study the possible formation of a medium of quarks and gluons under extreme 
temperature and energy density. Before the collisions, all the energy is carried by the beam particles in longitu- 
dinal phase space. But after the collisions, the final state particle production in the transverse plane carries finite 
transverse energy (Ey), which is an event-by-event observable and it is closely related to the collision geometry. 
In the Bjorken boost-invariant hydrodynamics model’ for relativistic heavy-ion collisions, Ey at mid-rapidity 
gives the quantitative estimation of the initial energy density produced in an interaction. Under boost invari- 
ance, the Bjorken energy density (gj) in the nuclear collision zone can be estimated as, 


1 dEr 


EBj = TS dy” (2) 


where, t is the formation time and usually taken to be 1 fm/c. Ey is the total transverse energy and Sy = 7 R? is 
the transverse overlap area of the colliding nuclei. As R = RoA?, one replaces A = Npart/2. That makes the 
expression for transverse overlap area, 
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Figure 2. Top (bottom) panel: transverse spherocity distribution for different centrality (multiplicity) classes in 
Pb-Pb (pp) collisions at ./snn = 5.02 TeV (./s = 13 TeV) using AMPT (PYTHIA8) model. 





Sy = 1h (Ae ) - (3) 
Transverse energy (Ey) at mid-rapidity region can be approximated as'?””', 
a dN dN 
ae aS (om) +2x (10) cao (4) 


The multiplicative factors 3/2 and 2 account for neutral particles. mp = ,/ Pe + m?, is the transverse mass 
and dN/dy is the integrated yield for 7*, K~ and p + p at mid-rapidity region i.e. |y| < 0.5, estimated for 
Pr >0.15 GeV/c. 

Figure 3 shows the mean transverse mass ((my)) as a function of different centrality classes in Pb-Pb collisions 
at ./SNN = 5.02 TeV at mid-rapidity for 7*, K* and p + pusing AMPT model. Except pions, where resonance 
decay contributions are expected, for all other charged particles (my) shows a decrease towards peripheral colli- 
sions. The lower panels of the figure show the effect of event topology on (my), where one observes higher (mr) 
for low-So events (jetty). Figure 4 shows the integrated yield (dN/dy) as a function of different centrality classes 
in Pb-Pb collisions at ,/synj = 5.02 TeV at mid-rapidity for different identified particles using AMPT model. As 
expected, the integrated yield is higher for more central collisions and gradually decreases as we move towards 
mid-central and peripheral collisions. It can be accounted due to the decrease in the participating partonic mat- 
ter from central to peripheral collisions. For pions, the integrated yield is higher than kaon and proton, which 
follows a thermalised Boltzmannian production of particles in a multiparticle production process. The integrated 
yield as a function of spherocity shows that high-So events have a higher yield than low-Sp events. It is also clear 
that the integrated yield is highly dependent on the spherocity classes for most central heavy-ion collisions and 
the dependence decreases while going towards peripheral collisions. 

Figure 5 shows the Bjorken energy density (¢g)) vs. centrality (%) for high-So, So-integrated and low-Sp events 
in Pb-Pb collisions. We observe a strong dependence of Bjorken energy density on the centrality classes. The 
values of the initial energy density is observed to be higher than the lattice QCD estimation of 1 GeV/fm? energy 
density for a deconfinement transition”. However, we have found that the Bjorken energy density is independent 
of the spherocity selection and irrespective of the event topology, it is similar for both high-Sp and low-Sy events 
for all collision centralities. It is noteworthy that the Bjorken energy density for a given collision centrality has 
effects from the number of particles produced and their mean transverse mass. Because of the opposite trends 
of these two observables with event spherocity, the overall effect cancels out leaving Bjorken energy density to 
be independent of event topology. 
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Figure 3. (mr) vs. collision centrality(%) for pions (top), kaons (middle) and protons (bottom) in Pb-Pb 
collisions with high-So, So-integrated and low-So events, respectively. 


Squared speed of sound (c2) and pseudorapidity distribution. Figure 6 shows pseudorapidity dis- 
tributions of charged particles in Pb-Pb collisions at ,/syn = 5.02 TeV at mid rapidity for (0-10)% and (60-70)% 
centrality classes in different spherocity classes. Figure 6 complements the observation seen in Fig. 4, where the 
charged particle multiplicity is found to be higher for high-Sp events compared to low-Sp events. 

In Landau hydrodynamical model”’, the width of the rapidity distribution is related to the speed of sound 


(cs) via the following expression. 
8 SNN 
a == a, ( YN), (5) 
31-—c2 2mp 
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Figure 4. Integrated yield (dN/dy) vs. collision centrality (%) for pions (top), kaons (middle) and protons 
(bottom) at mid-rapidity in Pb-Pb collisions. 
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Figure 5. Top plot: Bjorken energy density (¢g;) vs. centrality (%) with high-So, So-integrated and low-So 
events in Pb-Pb collisions. Bottom plot: ratio of Bjorken energy density (¢g;) for high-Sp and low-So events 


to the So-integrated events. The dotted line shows lattice QCD predicted value of critical energy density for a 
deconfinement transition. 


Here, mp is the mass of proton and oy is the width of the rapidity distribution and c? = 1/3 for ideal gas. Due 
to presence of a dip, it is difficult to fit a single Gaussian function to the pseudorapidity distribution. Generally, 
in experiments™ a double Gaussian function is used to fit the pseudorapidity distribution, which is given by, 


=1?2 =? 


Aye2 — Axe? A (6) 


Here, A; and A2 are normalisation parameters and o; and 02 are the widths of the double Gaussian distribution. 
After fitting we have obtained the values of oj and 02. The values of 0; and o2 are given in Table 2. The fitting 
has been performed using x” minimisation method and corresponding x*/ndf values for the fittings for each 
spherocity classes across all centralities are shown in Table 3. Here, o has been used as the default value and 
the maximum deviation of 7; and 02 is used for the uncertainty calculation of c?. Figure 7 shows the squared speed 
of sound as a function of centrality in Pb-Pb collisions at ./sjn = 5.02 TeV for different spherocity classes. The 
c* value as a function of centrality shows that for the central collision system, it is higher and gradually decreases 
towards peripheral collisions. It indicates that central collisions are denser compared to the peripheral collisions. 
However, c? is found to be similar for all spherocity classes within uncertainty. 


Kinetic freeze-out properties. The hot and dense medium formed in relativistic heavy-ion collisions 
cools down as the system expands until the kinetic freeze-out is achieved. At kinetic freeze-out, the transverse 
momentum spectra of the particles are frozen, which carries information about the phase-space distribution of 
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Figure 6. Charged particles pseudorapidity distribution (dN,;,/dn) for high-So, So-integrated and low-So events 
for (0-10)% (top) and (60-70)% (bottom) centrality class in Pb-Pb collisions at ./syjn = 5.02 TeV. 












































Centrality(%) | High-So So Integrated Low-So 
O71 02 O71 02 oO} 02 

0-10 2.110 + 0.024 | 1.757+0.223 | 2.120+0.023 | 1.786+0.022 | 2.079+0.022 | 1.769 + 0.020 
10-20 2.050 + 0.019 | 1.819+0.018 | 2.053+0.016 | 1.830+0.016 | 2.007+0.017 | 1.797 +0.016 
20-30 1.946 + 0.011 | 1.794+0.019 | 1.954+0.010 | 1.807+0.010 | 1.911+0.010 | 1.773 + 0.016 
30-40 1.845 + 0.006 | 1.751+0.006 | 1.853+0.006 | 1.759+0.006 | 1.830+0.006 | 1.740 + 0.006 
40-50 1.815 +0.004 | 1.747+0.004 | 1.803 +0.004 | 1.745+40.004 | 1.790+0.004 | 1.731 + 0.004 
50-60 1.755 +0.001 | 1.742+0.001 | 1.746+0.001 | 1.733+0.001 | 1.735+0.003 | 1.722 + 0.003 
60-70 1.738 + 0.034 | 1.731 40.034 | 1.743 40.003 | 1.737+0.003 | 1.729+0.003 | 1.722 + 0.003 
70-100 1.846 +0.001 | 1.843+0.001 | 1.799+0.001 | 1.798+0.001 | 1.782+0.001 | 1.781 +0.001 

















Table 2. Double Gaussian width parameters from fitting the pseudorapidity distributions in the range |n| < 2 
using Eq. (6). 


the final state of the fireball produced during relativistic heavy-ion collisions. The Boltzmann-Gibbs blastwave 
(BGBW) distribution” can be used to describe the transverse momentum spectra of identified charged particles 
formed in heavy-ion collisions and one can obtain the transverse radial flow velocity (fy) and kinetic freeze-out 
temperature (Tin) of the system. The invariant yield in the BGBW framework can be expressed as: 


a°N My 
EE, = c | ao.p"ep( - rn). (7) 


Here, C is the normalisation constant and p", the particle four momentum is given by, 
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x7/ndf 
Centrality(%) | High-So | So Integrated | Low-So 
0-10 0.16 0.11 0.13 
10-20 0.55 0.45 0.30 
20-30 0.92 0.99 1.15 
30-40 1.16 1.28 0.97 
40-50 1.06 LL7 0.92 
50-60 0.64 0.86 0.96 
60-70 0.67 0.64 0.91 
70-100 0.51 0.35 1.15 




















Table 3. x? /ndf values for the fitting of dN,;,/dn to a double Gaussian distribution. 
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Figure 7. Top plot: Speed of sound (c2) vs. Centrality(%) for high-So, So-integrated and low-Sp events in Pb—-Pb 
collisions. Bottom plot: Ratio of speed of sound (c2) for high-So and low-So events to the So-integrated events. 


p" = (mrcoshy, pr cos®¢, pr sing, mrsinhy). 
The particle four-velocity is given by, 
u“ = coshp (coshn, tanhp cos¢,, tanhp sing,, sinh 7), 
so that 
pup, = mr cosh(y — n) cosh p — pr sinh p cos(¢ — ¢). 
The freeze-out surface is parametrised as, 
Poy = (coshn, 0, 0,— sinh) t r dr dn d¢,. 


Here, 7 is the space-time rapidity. Now, Eq. (7) can be written as, 








aN Ro h inh 
P| epee fa xy (ME SMP) (PE she), 
dprdy y=0 0 Tkin Tin 
mr coshp pr sinhp . ; . ; . 
Here, K, (=) and Ip (a) are modified Bessel’s functions, which are given by, 
kin kin 
h oe h h 
Ky (— *) = | coshy exp( ae — » bash dy, 
kin 0 kin 
pr sinhp 1 [7 pr sinhp cosp 
(= —) = =f ex (=F jas, 
Tkin 2 Jo . Tkin 


(8) 


(9) 


(10) 


(11) 


(12) 
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Figure 8. Simultaneous fitting of identified charged particles’ py spectra with Boltzmann-Gibbs blastwave 
function for high-So (left), So-integrated (middle) and low-So (right) events in Pb-Pb collisions at ./snn = 5.02 
TeV. The fitting is shown for (0-10)% (top) and (60-70)% (bottom) centrality classes. 
































x7/ndf 
Centrality(%) | High-So | So Integrated | Low-So 
0-10 2.7 2.8 3:2. 
10-20 25 2.8 2.8 
20-30 2.9 2.8 2.7 
30-40 2.8 2.6 2.6 
40-50 2.4 2.2 2.1 
50-60 2.1 2.1 2.0 
60-70 2.4 22, 1.6 
70-100 2.2 2.1 2.4 




















Table 4. x7/ndf values for the simultaneous fitting of identified charged particles’ py spectra to the BGBW 
distribution. 


where p = tanh”! Brand Br = fsé™>*. Bris called radial flow velocity, € = (r/Ro), Bs is the maximum surface 
velocity, ris the radial distance and Ro is the maximum radius of the fireball at freeze-out. In this model, the 
particles closer to the center of the fireball are assumed to move slower than the ones at the edges. The mean 
transverse velocity is given by”, 


_ [Bere de 2 
(Or) = eae = (aya) Pe (13) 


Figure 8 shows the simultaneous BGBW fitting to the identified charged particles’ py spectra in (0-10)% and 
(60-70)% centrality classes. The fitting ranges for pions, kaons and protons are (0.5-1) GeV/c, (0.4-1.5) GeV/c 
and (0.3-3) GeV/c, respectively. The fitting has been performed using x” minimisation method keeping Tkin Bs 
and nas free parameters. Here, n designates the order of the flow profile. The corresponding x*/ndf values for 
the fittings for each spherocity classes across all centralities are shown in Table 4. 

Figure 9 shows the variation of kinetic freeze-out temperature (Tin) versus mean transverse radial flow 
velocity ((r)) for different spherocity and centrality classes obtained from simultaneous fit of identified charged 
particles’ py spectra with BGBW function for Pb-Pb collisions at ./syn = 5.02 TeV. The shaded area shows the 
uncertainties from the simultaneous BGBW fits. As we move from central to peripheral collisions, Tyin increases 
and (fr) decreases for all the spherocity classes. This behavior is naively expected due to the fact that central 
collisions are expected to have higher multiplicity in the final state, which would require more time to reach the 
freeze-out than the peripheral collisions. Conversely, with the increase in multiplicity one would expect higher 
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Figure 9. Kinetic freeze-out temperature versus mean transverse radial flow velocity obtained from 
simultaneous fit of identified particles’ py-spectra with BGBW distribution function for high-So (blue circle) 
and low-So (red triangle) in Pb-Pb collisions at ./snn = 5.02 TeV. Linear fits for high-So and low-So are shown in 


dotted blue and red lines, respectively. The shaded area shows the uncertainties from the simultaneous BGBW 
fits. 
































High-So So Integrated Low-So 

Centrality(%) _| (Br) Tkin [GeV] _| (Br) Tkin[GeV] | (Br) Tkin [GeV] 

0-10 0.539 + 0.007 | 0.124+0.002 | 0.543+0.006 | 0.123+0.002 | 0.545+40.007 | 0.124+0.002 
10-20 0.550 + 0.007 | 0.123+0.002 | 0.549+0.006 | 0.123+0.002 | 0.540+0.008 | 0.129 + 0.003 
20-30 0.549 + 0.006 | 0.122+0.002 | 0.553+0.006 | 0.122+0.002 | 0.549+0.007 | 0.128 + 0.003 
30-40 0.548 + 0.007 | 0.120+0.002 | 0.546+0.006 | 0.123+0.002 | 0.542 +0.007 | 0.128 + 0.003 
40-50 0.538 + 0.006 | 0.121 +0.002 | 0.535+0.006 | 0.125+0.002 | 0.532+0.007 | 0.130 + 0.002 
50-60 0.519+0.007 | 0.125+0.002 | 0.515+0.007 | 0.129+0.002 | 0.514+0.008 | 0.134 + 0.003 
60-70 0.492 + 0.008 | 0.129+0.003 | 0.484+0.008 | 0.134+0.003 | 0.483+0.009 | 0.139 + 0.003 
70-100 0.437 40.010 | 0.136+0.003 | 0.407+0.015 | 0.145+0.007 | 0.307+0.021 | 0.163 + 0.008 





























Table 5. Kinetic freeze-out temperature (Tj) and mean transverse radial flow velocity ((6r)) obtained from 
simultaneous fit of identified charged particles’ py-spectra with Boltzmann-Gibbs blastwave function. 


transverse radial flow, which we observe in Fig. 9. The spherocity dependent values of (By) and Tin for all cen- 
tralities for Pb-Pb collisions at ./snn = 5.02 TeV are enlisted in Table 5. 

Contrary to the observables studied so far in the previous sections, we see a clear dependence of the kinetic 
freeze-out parameters on So. An isotropic (high-So) event is expected to be dominated by large number of soft 
particles, which would require more time to reach the freeze-out. Thus, the isotropic events are found to have 
lower Tyin compared to low-So events. To understand the dependence of the kinetic freeze-out parameters, we 
have fitted the parameters with a first-order polynomial function separately for each spherocity class. As evident 
from Fig. 9, the role of spherocity plays a bigger role in peripheral events and it diminishes when one goes towards 
central collisions. This could be an indication that spherocity is crucial while studying the final state effects. 


Summary and conclusion 

We have implemented transverse spherocity in Pb-Pb collisions at ./syn = 5.02 TeV using AMPT model and 
study the dependence of transverse spherocity on different global observables in heavy-ion collisions at the 
Large Hadron Collider energies. In summary, the spherocity distributions in Pb-Pb collisions are found to be 
shifted more towards the isotropic limit when compared to the pp collisions, where the distributions are shifted 
towards the jetty limit. This behavior is understood based on the fact that the system size in Pb-Pb collisions are 
significantly higher when compared with pp collisions and because of the medium effects through the process of 
isotropisation, the jettiness of the events gets suppressed to a larger extent. The Bjorken energy density and speed 
of sound are found to be independent of the spherocity selection in heavy-ion collisions. However, we found 
that kinetic freeze-out parameters depend on spherocity. The role of spherocity plays a bigger role in peripheral 
events and it diminishes when one goes towards central collisions. The sensitivity of event topology however, 
depends on the observable under study because of some of the counter-balancing effects in view of medium 
effects in heavy-ion collisions or high-multiplicity environments. 





Scientific Reports | 


(2022) 12:3917 | https://doi.org/10.1038/s41598-022-07547-z nature portfolio 


www.nature.com/scientificreports/ 


Received: 8 December 2021; Accepted: 15 February 2022 
Published online: 10 March 2022 


References 


1; 


2. 


3; 


4, 


5, 


Mallick, N., Sahoo, R., Tripathy, S. & Ortiz, A. Study of transverse spherocity and Azimuthal anisotropy in Pb-Pb collisions at 
./SNN = 5.02 TeV using A Multi-Phase Transport Model. J. Phys. G 48, 045104. https://doi.org/10.1088/1361-6471/abeb59 (2021). 
Lin, Z.-W., Ko, C. M., Li, B.-A., Zhang, B. & Pal, S. A Multi-phase transport model for relativistic heavy ion collisions. Phys. Rev. 
C72, 064901. https://doi.org/10.1103/PhysRevC.72.064901 (2005). 

Wang, X.-N. & Gyulassy, M. HIJING: A Monte Carlo model for multiple jet production in p p, p A and A A collisions. Phys. Rev. 
D 44, 3501. https://doi.org/10.1103/PhysRevD.44.3501 (1991). 

Zhang, B. ZPC 1.0.1: A Parton cascade for ultrarelativistic heavy ion collisions. Comput. Phys. Commun. 109, 193. https://doi.org/ 
10.1016/S0010-4655(98)00010-1 (1998). 

Greco, V., Ko, C. M. & Levai, P. Parton coalescence and anti-proton/pion anomaly at RHIC. Phys. Rev. Lett. 90, 202302. https:// 
doi.org/10.1103/PhysRevLett.90.202302 (2003). 


. Greco, V., Ko, C. M. & Levai, P. Parton coalescence at RHIC. Phys. Rev. C 68, 034904. https://doi.org/10.1103/PhysRevC.68.034904 


(2003). 


. Li, B., Sustich, A. T., Zhang, B. & Ko, C. M. Studies of superdense hadronic matter in a relativistic transport model. Int. J. Mod. 


Phys. E 10, 267. https://doi.org/10.1142/S0218301301000575 (2001). 


. Li, B.-A. & Ko, C. M. Formation of superdense hadronic matter in high-energy heavy ion collisions. Phys. Rev. C 52, 2037. https:// 


doi.org/10.1103/PhysRevC.52.2037 (1995). 


. Fries, R. J., Miiller, B., Nonaka, C. & Bass, S. A. Hadronization in heavy-ion collisions: Recombination and fragmentation of partons. 


Phys. Rev. Lett. 90, 202303. https://doi.org/10.1103/PhysRevLett.90.202303 (2003). 


. Fries, R. J., Muller, B., Nonaka, C. & Bass, S. A. Hadron production in heavy ion collisions: Fragmentation and recombination 


from a dense parton phase. Phys. Rev. C 68, 044902. https://doi.org/10.1103/PhysRevC.68.044902 (2003). 


. Tripathy, S., De, S., Younus, M. & Sahoo, R. Predictions for azimuthal anisotropy in Xe+Xe collisions at ,/syn = 5.44 TeV using 


a multiphase transport model. Phys. Rev. C 98, 064904. https://doi.org/10.1103/PhysRevC.98.064904 (2018). 


. Loizides, C., Kamin, J. & d’Enterria, D. Improved Monte Carlo Glauber predictions at present and future nuclear colliders. Phys. 


Rev. C 97, 054910. https://doi.org/10.1103/PhysRevC.97.054910 (2018). [Erratum: Phys. Rev. C 99, 019901 (2019)]. 


. Cuautle, E. et al. Disentangling the soft and hard components of the pp collisions using the sphero(i)city approach (2014). arXiv: 


1404.2372 [hep-ph . 


. Ortiz, A., Paic, G. & Cuautle, E. Mid-rapidity charged hadron transverse spherocity in pp collisions simulated with Pythia. Nucl. 


Phys. A 941, 78. https://doi.org/10.1016/j.nuclphysa.2015.05.010 (2015). 


. Salam, G. P. Towards Jetography. Eur. Phys. J. C 67, 637. https://doi.org/10.1140/epjc/s10052-010-1314-6 (2010). 
. Sjéstrand, T. et al. An introduction to PYTHIA 8.2. Comput. Phys. Commun. 191, 159. https://doi.org/10.1016/j.cpc.2015.01.024 


(2015). 


. Khuntia, A., Tripathy, S., Bisht, A. & Sahoo, R. Event shape engineering and multiplicity dependent study of identified particle 


production in proton + proton collisions at ./s = 13 TeV using PYTHIA8. J. Phys. G 48, 035102. https://doi.org/10.1088/1361-6471/ 
abb1f8 (2021). 


. Bjorken, J. D. Highly relativistic nucleus-nucleus collisions: The central rapidity region. Phys. Rev. D 27, 140. https://doi.org/10. 


1103/PhysRevD.27.140 (1983). 


. Adam, J. et al. Measurement of transverse energy at midrapidity in Pb-Pb collisions at ./syn = 2.76 TeV. Phys. Rev. C 94, 034903. 


https://doi.org/10.1103/PhysRevC.94.034903 (2016). 


. Sahoo, R., Mishra, A. N., Behera, N. K. & Nandi, B. K. Charged particle, photon multiplicity, and transverse energy production in 


high-energy heavy-ion collisions. Adv. High Energy Phys. 2015, 612390. https://doi.org/10.1155/2015/612390 (2015). [Erratum: 
Adv. High Energy Phys. 2017, 4517153 (2017)]. 


. Abelev, B. I. et al. Systematic measurements of identified particle spectra in pp, d* Au and Au+Au collisions from STAR. Phys. 


Rev. C79, 034909. https://doi.org/10.1103/PhysRevC.79.034909 (2009). 


. Karsch, F. Lattice results on QCD thermodynamics. Nucl. Phys. A 698, 199. https://doi.org/10.1016/S0375-9474(01)01365-3 (2002). 
. Landau, L. D. On the multiparticle production in high-energy collisions. Izv. Akad. Nauk Ser. Fiz. 17, 51 (1953). 
. Abbas, E. et al. Centrality dependence of the pseudorapidity density distribution for charged particles in Pb-Pb collisions at ./sNn 


= 2.76 TeV. Phys. Lett. B 726, 610. https://doi.org/10.1016/j.physletb.2013.09.022 (2013). 


. Schnedermann, E., Sollfrank, J. & Heinz, U. W. Thermal phenomenology of hadrons from 200-A/GeV S+S collisions. Phys. Rev. 


C 48, 2462. https://doi.org/10.1103/PhysRevC.48.2462 (1993). 


. Huovinen, P., Kolb, P. E, Heinz, U. W., Ruuskanen, P. V. & Voloshin, S. A. Radial and elliptic flow at RHIC: Further predictions. 


Phys. Lett. B 503, 58. https://doi.org/10.1016/S0370-2693(01)00219-2 (2001). 


. Braun-Munzinger, P., Stachel, J., Wessels, J. P. & Xu, N. Thermal equilibration and expansion in nucleus-nucleus collisions at the 


AGS. Phys. Lett. B 344, 43. https://doi.org/10.1016/0370-2693(94)01534-J (1995). 


. Tang, Z. et al. Statistical origin of constituent-quark scaling in the QGP hadronization. Chin. Phys. Lett. 30, 031201. https://doi. 


org/10.1088/0256-307X/30/3/031201 (2013). 


. Adcox, K. et al. Single identified hadron spectra from s(NN)**(1/2) = 130-GeV Au+Au collisions. Phys. Rev. C 69, 024904. https:// 


doi.org/10.1103/PhysRevC.69.024904 (2004). 


Acknowledgements 

RS acknowledges the financial supports under the CERN Scientific Associateship and the financial grants under 
DAE-BRNS Project No. 58/14/29/2019-BRNS of Government of India. DB acknowledges the financial supports 
from CSIR, Government of India. SP acknowledges the financial supports from UGC, Government of India. ST 
acknowledges the supports under the INEN postdoctoral fellowship. The authors would like to acknowledge 
the usage of resources of the LHC grid computing facility at VECC, Kolkata and usage of resources of the LHC 
grid Tier-3 computing facility at IIT Indore. The authors acknowledge Dr. Arvind Khuntia and Mr. Dushmanta 
Sahu for careful reading of the manuscript. 


Author contributions 
All authors have equally contributed to the research article in terms of formulating the problem through discus- 
sions, data generation, data analysis, manuscript preparation, editing and correspondence. 


Competing interests 
The authors declare no competing interests. 





Scientific Reports | 


(2022) 12:3917 | https://doi.org/10.1038/s41598-022-07547-z nature portfolio 


www.nature.com/scientificreports/ 


Additional information 
Correspondence and requests for materials should be addressed to R.S. 


Reprints and permissions information is available at www.nature.com/reprints. 


Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and 
institutional affiliations. 


Open Access This article is licensed under a Creative Commons Attribution 4.0 International 
BY License, which permits use, sharing, adaptation, distribution and reproduction in any medium or 


format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the 
Creative Commons licence, and indicate if changes were made. The images or other third party material in this 
article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the 
material. If material is not included in the articles Creative Commons licence and your intended use is not 
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from 
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 


© The Author(s) 2022 





Scientific Reports | (2022) 12:3917 | https://doi.org/10.1038/s41598-022-07547-z nature portfolio 


